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SUMMARY 


cin attempt is made to develop a second approximation to the 
solution of problems of supersonic flew which can be solved by 
existing first-order theory. The method of attack adopted is an 
iteration process using the linearized solution as the first step. 

For plane flow it is found thai a particular integral of the 
iteration equation can be writien down at once in terms of the 
first-order solution. The second-order problem ta thereby 
reduced to an equivalent first-order problem and can be readily 
solved, At the surface of an isolated body, the solution reduces 
io the well-known result of Busemann. The plane case ie 
considered in some detail insofar as tt gires insight into the 
nature of the iteration process. 

Again, for anally symmetrie flow the problem is reduced to a 
first-order problem by the discorery of a particular integral. 
For smooth bodies, the second-order solution can then be cal- 
culated by the method of ron Kármán and Moore. Bodies 
with corners are also treated by a slight modification of the 
method. The second-order solution for cones represents a 
considerable tmprorement over the linearized result. Second- 
arder theory also agrees well with seceral solutions for other 
bodies of revolution calculated by the numerical method of 
characteristics. 

For full three-dimensional flow, only a partial particular 
integral has been found. As an example of a more general 
problem, the solution is dericed for an inclined cone. The 
possibility of treating other inclined bodies of revolution and 
three-dimenstonal wings is discussed briefly. 


INTRODUCTION 


As the linearized theory of supersonic flow approaches full 
development, the question arises whether more exact approx- 
imations are practical. If viscous effects are large, refine- 
ment of the perfect-fluid solution is impractical. If viscosity 
is negligible, however, higher approximations are known to 
yield a closer approach to reality. In intermediate cases, 
un improved solution is desirable in order to assess the 
relative effects of viscosity and nonlinearity. 

The prototype of a higher-order solution for supersonic 
fow is Busemann’s series for the surface pressure in plane 
flow past an isolated body. This simple result is of con- 
siderable value in analyzing supersonic airfoil sections. 
Two terms of the series prove sufficient for almost all require- 
ments; the extension to third and fourth order is chiefly of 
academic interest. 

The aim of the present study is, therefore, to find a second 
approximation, analogous to Busemann’s result, for super- 


sonic flow past bodies which can be treated by existing 
first-order theory. The natural method of attack, and 
apparently the only practical one, 1s by means of an iteration 
process, taking the usual linearized result as the first step. 
Several writers have applied this procedure to two-dimen- 
sional subsonic flow. In supersonic flow, as usual, the 
solution is simpler, so that more general problems can be 
solved. 

This paper is based upon a thesis for the degree of doctor 
of philosophy in eeronautics written at the California 
Institute of Technology under a National Research Council 
predoctoral fellowship and under the guidance of Prof. P. A. 
Lagerstrom (reference 1). It was published m revised form 
as NACA TN 2200, 1951 (reference 2). The present version 
has been further slightly revised, in particular to include ref- 
erences to the recent literature. 


ITERATION PROCEDURE 
BASIC ASSUMPTIONS 


The problem to be considered is that of steady three- 
dimensional supersonic flow of a polytropic gas past one or 
more slender bodies. As indicated in figure 1, the bodies 





Fiter {.—The problem. 


are assumed either to be pointed or to extend upstream 
indefinitely as cylinders parallel to the free-stream direction. 
In either case, the origin of coordinates can be chosen so 
thet all variations in body shape are confined to the half- 
space r>>0. Wind axes are introduced, so that for z<0 the 
flow is uniform and parallel to the z axis, with the freestream 
velocity U and Mach number M. (For definitions of all 
symbols, see the appendix. ) 

The bodies are slender, which means that at any point the 
component of LU normal to the surface is small compared 


1Supersedes NACA TN 200, “4 Study of Second-Order Supersonic-Flow Theory,” by Milton D. Yan Dyke, 1951. 
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with U itself. The symbol e will be used throughout as a 
measure of this smallness. Thus the ordinates of a body 
will be written as e times a function of order unity. Used 
in this way, e serves to distinguish terms of various orders. 

The aim of this investigation is to find a second approxima- 
tion to the solution of problems for which the first-order 
solution is available. The first-order, or linearized, solution 
is defined as the result of keeping only linear perturbation 
terms in the equation of motion. 
solution is the result of retaining second-degree terms .in 
perturbation quantities. In addition, however, certain of the 
triple products are in some cases found to be as important 
as one or more double products and are therefore also re- 
tained in the equation. It may be noted that the second- 
order solution will not generally consist simply of terms of 
order e and è, though this is the case for plane flow. For 
example, the second-order solution for flow past a body of 
revolution will be found to contain terms as high as éine. 

The flow is assumed to be irrotational and isentropic. 
This assumption is justified in the first- and second-order 
solutions, since the resulting error is found to be at most of the 
order of terms neglected elsewhere. 


EXACT PERTURBATION EQUATION 
Under the assumption of irrotational flow, there exists a 
velocity potential Q. In Cartesian coordinates, the equation 
of motion is (reference 3, equation (39)). 
(e —27)922+ (?’—0,7)Q,,+ (ce? —0,7)0,,— 
22,.2,2yp— 20,0,.022— 20,0, 7)==0 (1) 


Here the local speed of sound ¢ is related to c., its value in 
the uniform stream, by 





Pme 1 (0,7+-2,7+-0,?— U’) (2) 


where y is the adiabatic exponent. The subscript notation 
is used to indicate differentiation. 

A perturbation potential # is now introduced in the usual 
way. For convenience, however, ® is normalized through 


division by the free-stream velocity so that 
a= U(r+) (3) 


The perturbation velocity at any point is then the gradient 
of multiplied by U. - 

Introducing the perturbation potential into the equation 
of motion gives, after some manipulation, 


byt Pu BBM | T+ 28,4824 OF + 82) Oat 


Dy + Buz) +28, Fee + Dat nt 
b26,,+28,8,8,,+28,(14+6,)8,.+ 


2(1 +8.) Obey | (4) 


where @=/A?—1. 
SOLUTION BY ITERATION 


The exact perturbation equation (equation (4)) 1s com- 
pletely equivalent to the original nonlinear potential equa- 
tion (equation (1)). Simplifying assumptions must there- 


Similarly, the second-order . 
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fore be introduced in order to solve it. If it is assumed 
that squares and products of the derivatives of @ can be 
neglected, the right-hand side of equation (4) disappears, 
leaving the wave equation l 


6,0 +,,°) me Bd, =0 (5) 


This equation is the basis of the linearized or first-order 
perturbation theory, and its solution is designated by 6”, 

Nore exact solution of equation (4) by means of iteration 
was first suggested by Prandtl (reference 4). The procedure 
has been applied to plane subsonic flow by Gértler (refer- 
ence 5), Hantsche and Wendt (references 6 and 7), Imai and 
Oyama (references 8 to 10), and Kaplan (references 11 to 
13}. Schmieden and Kawalki (reference 14) applied the pro- 
cedure to subsonic flow past an ellipsoid of revolution. Most 
of these writers have considered the stream function rather 
than the potential, which restricts the method to plane or 
axially symmetric flows. The procedure is clearly described 
by Sauer (reference 3, p. 140) for the case of plane flow. 

The linearized solution 6%, subject to proper boundary 
conditions, is taken as the first approximation. Substi- 
tuting this known solution into the right-hand side of 
equation (4) gives 


Pr” tOu” — pbu =F; (2,y,2) (6) 


where F is a known function of the independent variables. 
This is again a linear equation, the nonhomogeneous wave 
equation. A second-order solution ®, subject to. proper 
boundary conditions, can be sought by standard methods. 
The procedure can be repeated by substituting 6™ into the 
right-hand side of equation (4) and solving again. Con- 
tinuing this process yields a sequence of solutions J% 
which, under proper conditions, perhaps converges to the 
exact solution. 

A significant feature of this procedure is thet in each step 
the left-hand side of the iteration equation js the same. As 
a consequence, the characteristic curves of each iteration 
equation are just the Mach lines of the undisturbed flow. 
However, in actuality the local Mach lines are usually 
neither straight nor parallel; that is, the characteristics of 
the original nonlinear equation are curved, in a manner which 
is initially unknown and which depends upon the shape of 
the body. Because of the fundamental role played by the 
characteristics in the theory of hyperbolic equations (see, 
e. g., reference 15, ch. 5; reference 16, ch. 2), it might be 
anticipated that an iteration procedure should be chosen 
such that in each step the approximate characteristics 
would be successively revised so as to approach the actual 
characteristics. For purely subsonic flow, the counterpart 
of such & procedure is known to converge under proper 
conditions (reference 15, p. 288-289). Convergence might 
reasonably be anticipated also in the case of purely super- 
sonic flow. 

Unfortunately, an iteration procedure in which the 
approximate characteristics are successively revised would, 
except in the first step, involve equations with nonconstant 
coefficients. This would greatly complicate the procedure. 
Fortunately, it will be found that the scheme adopted here, 
which makes no provision for such revision, nevertheless 
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gives an improved solution nearly everywhere in the flow 
feld. 


SECOND-ORDER ITERATION EQUATION 


Henceforth, only the first two steps of the iteration process 
will be considered in detail. In order to eliminate cumber- 
some superscripts, it is therefore convenient to introduce & 
simplified notation for the first two approximations. The 
first-order perturbation potential will hereafter be denoted 
by g, and the second-order potential? by ¢ 


g=p ) 


m (7) 


Introducing these quantities into the exact perturbation 


equation (equation (4)) gives the following second-order 
iteration equation for ¢: 





Pyr t $u Bpr =A |? = E (2e.+ ee test gs") (gzz + a t 


Puzl +2 PrPrr + Pr Ors + Pr Spy tT Pr Pz + 
Deffa Fell + G2) Giz + 2(14+ e) oron | (8) 


Since ¢ satisfies equation (5), the term (pret purt fz) 0 
the right-hand side of equation (8) can be replaced by 
Pes, and the equation for ¢ becomes 


Les 
5 AP Qe,+e7%+of7t+ pr) Prr T 





Dyy le Èz — B’prr = E 


2 Crrt Os Prt Pr Lurt Os Pet 
Ipu F 2¢.(1 +92) G2 +2(1 +2) Prez | 


(9) 


Here the right-hand side contains not only double products 
of perturbation quantities but also triple products. The 
latter can be omitted for plane flow, since they contribute 
terms of smaller order (equal to those found in the next 
iteration}. Otherwise, certain of the triple products should 
be retained because their contribution is as great as that 
of one or more of the double products and greater than any 
contribution from a third approximation. It will be seen 
later that triple products should be retained if they involve 
only derivatives normal to the free stream. Triple products 
which involve x derivitives can be neglected, so that the 
iteration equation becomes 


Prat pa — Pp =M [2 +(y— DAP] orp +2 eyezy T 
2 Paver + Py Pyy +2 PrP t Pr G22 j (10 ) 


Here the triple products which may be important are the 
last three terms on the right-hand size. 


The adiabatic exponent y will be found to occur always in 
the combination 
yar 1)AL° 


ae ya 


(11) 





t Here ġ is regarded as being the complete second-order potential rather than (as fn ref- 
erences 1 and 2) a small correction to be added to the first-order solution. 
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Introducing this expression in place of y gives the final form 
of the second-order iteration equation: 


rr t bu— Pbr =IP(20V— 1) Pererr 2P Fent 
Or Cys +2 g,020y2+ Pr Prz (12) 
TFERATION EQUATION IN OTHER COORDINATES 


In cylindrical coordinates, equation (9) becomes 


bt EHIS ppu IP] AN- DPen t Bee +2 E pat 








2 2 
Orpi Greve re 
Or Orn t2 2 Prt— i + rt geet 


of v2 2 g PrF 
Pr Prr! Pr Prr! ri Prr! r3 Cré? PrP Sarr 
(13) 


The terms whose form is indicated at the end of the equation 

are the triple products which will be found to be negligible. 

For conical flows it is convenient to introduce nonorthogo- 

nal conical coordinates (x, t, @) where 
Êr 

=f (14) 

If the body itself is conical, the perturbation potentials are 


reduced to functions of two variables by introducing conical 
perturbation potentials (reference 17} so that 


plz, i, 0 =x, 0) : (15) 


with a corresponding definition for $. The derivatives are 
given by 


APEJE ASA P = 
r= — tpr P= Pi r= ppa 
7 B an ee 
pr= Pp: Pr Pit Cor = pie — ipt (16) 
= _ BE 
pe=Tr se pa =T pet Irr = aa Pit 


with the same relations connecting ¢ and ¢. The iteration 
equation (equation (13}) becomes 


(1-8) put HEM N-DE) Pu 2B But 


z 
E gı lgıe— tpt 8 p put 2 5 PpD — 


p? — bh — 4 — — 
a pipe tog pront Oj Elie it 
= Se eet ee ce = 
Poo r:, Geert; B (g— tpi) pe (ce—tpe), 
AE A (17) 


Here the grouping of terms corresponds to that in equa- 
tion (13). 
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BOUNDARY CONDITIONS 


Physical considerations suggest that the flow should satisfy 
the following conditions: 

Tangency condition: The resultant flow is tangent to the 
surface of the body. 

Upstream condition: All flow perturbations vanish every- 
where upstream of the body. 

These two requirements are sufficient to determine the 
solution. The first imposes one condition along the time- 
like surface of the body, and the second may be regarded as 
imposing two conditions on a spacelike surface. This is the 
case of mixed boundary conditions (reference 15, p. 172) 
and leads to a determinate solution (see reference 16, p. 85). 

The tangency condition may be written formally as 

grad Q-grad S=0 at S=0 (18) 


where S=0 is the equation of the surface of the body. In 
a more useful form it becomes, for the first- and second-order 
problems, 


Pe = (slope) (1 + Pz) 
e= (slope) (1+ ¢:) 


(on the surface) (19a) 


(on the surface) (19b) 
Here ¢ means the cross-wind component of the normal 
derivative of ¢ at the surface of the body. In plane flow 
$: is $,, and in axially symmetrie flow ¢, is ¢,. The slope of 
the body is measured with respect to the free-stream 
direction. 
In first-order theory, the tangency condition is frequently 
approximated by neglecting ¢, compared with 1 in equation 
(19a), which causes only a second-order error. Correspond- 
ingly, in the second-order tangency condition (equation 
(19b)), ¢- can be replaced by its first-order counterpart gs 
with only a third-order error... Thus the tangency condition 
simplifies to 


(20a) 
(20b) 


v.=slope (on the surface) 


e= (slope). (1+ ¢-) 


This approximation will not be made except for plane flow. 
Another approximation in the tangency condition can be 


made for planer bodies. A planar body is one whose entire 
surface lies near to @ plane parallel to the free stream, say 
the plane y=0 (reference 17, p. 52). Thin flat wings are 
planar bodies, whereas slender pointed bodies of revolution 
are not. For a planar body the first-order tangency condi- 
tion can be imposed at the plane rather than on the surface 
of the body. Correspondingly, the second-order tangency 
condition can also be imposed at the plane provided that the 
difference in ¢, between the plane and the surface of the 
body is accounted for approximately by retaining the second 
term in its Taylor series expansion about y=0. Further- 
more, to first order g. is ø, for a planar body. Thus if the 
surface of a planar body is given by y= Y(z,2), the simplified 
tangency condition Is 


ey F,(1 zg Pr) 
+ Vow— Yip: = 


(on the surface) 


at y=0 (21a) 


Y,(1+¢-) at y=0 (21b) 
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The term. containing ¢, in equation (21b) accounts for the 
fact that to second order the cross-wind component ġe may 
not be vertical. Corresponding simplifications can be made 
for any quasi-cylindrical body, which is a body whose entire 
surface lies near to.a cylinder parallel to the free stream. 

Finally, the approximation which led to equations (20) 
may be adopted in addition to that just discussed for planar 
systems, in which case the tangency condition simplifies 
further to 


Py Je l 
p= Y,;( 1+ Pr) —Y py + Yip: 


The upstream condition implies that both ® and ®, vanish 
at the plane z=0, which completes the boundary conditions 
required. For the first- and second-order problems the 
upstream condition is therefore 


at y=0 (22a) 


at y=0 (22b) 


at z=0 (23a) 


(23D) 


p= pz= Ù 
p= h: = at r=0 


DETERMINATION OF PRESSURE 


When the potential field has been determined, the net 
velocity g at any point is given by 


.  F=(U+u)*+e?4+ (24) 
where 
u 
gT? 
v " 
U le, (25) 
®, 
W 
U + &, 


in Cartesian and cylindrical coordinates, respectively. Be- 
cause the flow is assumed to be isentropic, the pressure 
coefficient is given by 


pe are 


where p, and p, are the free-stream pressure and density. 

It is the practice in linearized theory also to simplify the 
pressure relation. Substituting equation (24) into equation 
(26) and expanding gives 


C ee 





1} (26) 











Wn gt v a uv rwt Af? aa 
C,=—2 Di +8, wat Mee tk ang a Ga ee 
3 2 2 2 p? 
EEEE wP) ee 


All the terms shown here explicitly may give contributions of 
second order; the remaining terms whose. form is merely 
indicated are always of higher order. In linearized theory 
only the first term is ordinarily retaincd. This is satisfactory 
for plane flow or flow past planar systems, since the contribu- 
tion of the remaining terms is definitely of higher order. 
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u v? 


————— First order, Cp =-2>--= 
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Exact (reference 20) 


= +B | m 
-i z > 
Lán ut(1--#) T 3) 








FisURe 2.— Effect of using various pressure relations in first-order sotution for 5° cone. 


In fact, for plane flow past an isolated body it happens that 
the next two terms cancel identically. 

For slender bodies such as a cone, however, orders of 
magnitude are not so clearly distinguished. Busemann 
suggests (reference 18) that the second term, (+t w/t, 
is then sufficiently large compared with the first that it 
should be retained also. This view is supported by Lighthill 
(reference 19), who shows that the resulting solution is 
correct up to the order of the quantities contributed by the 
second term. Again, the third term, 67u?/L", which is also 
the square of a perturbation quantity, is comparable with 
the second at high Mach numbers (where Lighthill's order 
estimates become invalid) and might logically be retained. 
Having gone this far, it may be simpler to use the exact 
relation. | 

Each of these four possibilities for the first-order flow past 
à 5° cone is compared with the exact solution (reference 20) 
in figure 2. The series (equation (27)) is seen to alternate in 
this case. It converges so slowly, however, that at moderate 
Mach numbers (say, near M=42) where first-order theory 
is most accurate, linearizing the pressure relation introduces 
much greater errors than linearizing only the equation of 
motion. Adding each of the quadratic terms in turn causes 
changes nearly as great as the error due directly to non- 
linearity in the equation. 

The point of view to be adopted here is that calculating 
the velocities and calculating the pressure are two essentially 


distinct operations. <A certain degree of approximation may 
be necessary in order to solve for the velocities, but the 
pressure relation need not then be approximated to the same 
extent simply for the sake of consistency. Forit may happen 
that the resulting errors (though of the same mathematical 
order) sre greater than those due to the original approxima- 
tion. Indeed, this is evidently the case at moderate Mach 
numbers in the first-order solution for a cone and will be 
found true to a greater extent in the second-order solution. 
Moreover, in the second-order solution so many terms of 
equation (27) must be retained that it is usually simpler to 
use the exact relation. For these reasons, the exact pressure 
equation (equation (26)) will be used throughout except in 
the case of plane flow. 
ROLE OF A PARTICULAR INTEGRAL 

The second-order iteration equation can be attacked by 
standard methods, and in the case of plane flow a solution 
can be found directly. For plane and axially symmetric 
flows, however, tt will be found that a particular integral of 
the iteration equation can be written down at once in terms 
of the first-order solution. This solves the problem because 
the complete solution consists of a particular integral plus a 
solution of the homogeneous equation, and the latter can be 
obtained by existing methods. That is, the second-order 
potential may be written as | 


o=v+x (28) 


where 
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y is any particular integral of the nonhomogeneous iteration 
equation 

x is the complementary function which serves to correct the 
tangency condition 


The problem for x is the usual first-order problem for which 
methods of solution are assumed to be known. 

The role of the particular integral is to transfer the non- 
homogeneity in the problem from the equation, where it is 
troublesome, to the boundary conditions, where it can be 
handled by existing theory. For linear partial differential 
equations it is always possible in principle to transfer non- 
homogencities in this way from the equation to the boundary 
conditions, and vice versa, by adding a suitable function to 
the dependent variable (see reference 21, p. 236). 

Since the particular solution ¥ will be found in terms of the 
first-order solution, it will vanish upstream of the plane z=0. 
Then the complementary function must also vanish there, so 
that the upstream condition for x 1s 


x(0, ¥, 2) =x2(0, y, 2) =0 (29) 


From equation (19b}, the tangency condition for x is found 


to be 
Pet xe= (slope) (1 Hyr tX) 


or, in the case of a planar body given by y= F(z, z), from equa- 
tion (22b), 


pt X= 71+ e) —Y prt Ysy: at y=0 (31) 


It should be noted that, although ¢ is of the same magni- 
tude as g, this is not necessarily true of either y or x alone. 


PLANE FLOW 


The second-order solution for conditions at an isolated sur- 
face in plane supersonic flow was given by Busemann (refer- 
ences 22 and 23). By using the iteration procedure, the 
solution will now be found throughout the flow field, includ- 
ing the case when several bodies interact. 

The solution for plane flow is of interest chiefly insofar as 
it serves as & guide in more complicated problems. In par- 
ticular, it provides insight into such details of the iteration 
process as the question of its success and the effect of sharp 
corners. 


(on the surface) (30) 


PARTICULAR INTEGRAL FOR PLANE FLOW 
The first-order equation for plane flow is 


Prr — B pr = 0 (32) 
The general solution is 
o(2,y) =h(a— By) +3 (et By) (33) 
where Å and 7 are functions determined by the first-order 
boundary conditions. i 


In the iteration equation, all triple products can be neg- 
lected, and equation (12) becomes 


Pyy— B'h =2M? [(V— 1) B Psp t PuPzy] (34) 


It can be verified directly that a particular integral of this 
equation is given in terms of the first-order solution by- 


=Ma | (1-3) e+ 5 ve | (35) 
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To this must be added & solution x of the homogeneous 
equation (equation (32)), which has the form 


x=H(z— By) + S(t By) (36) 


where H and J are functions determined by the second-order 
boundary condition. 

For flow past a single boundary (such as one surface of an 
airfoil) the first-order potential (equation (83)) contains 
only one or the other of the functions A and 7. In this case 
PyCry= Boe, SO that the iteration equation reduces to 


Pyy— B'a = 2AL RN Px Prs (3 7) 


The particular integral may then be simplified to 


N 
and the complementary function contains only H or J, 
according as the first-order solution contains only A or 7. 


FLOW PAST A CURVED WALL 


As an example of the application of the particular solution, 
consider fow past a wall which at some point begins to 
deviate slightly from a plane (see fig. 3). The wall can be 
represented by 


y= Y (x) =¢G(2) (39) 


where e is a parameter small compared with unity and g(x) 
is a continuous function of order unity which vanishes for 
sso. , : 
This is a planar body, so that the tangeney condition is 
given by equation (22a). Consequently, the first-order 
problem is 


Cy — Bos 
y(t, 0)=eg"(z) (40) 
90, y)=¢.(0, y) =0 
where the prime indicates differentiation of g with respect to 
its argument. The solution is 


p=—5 g(z— By) (41) 


Substituting this first-order solution into the right-hand 
side of equation (37) gives the iteration equation 


by — Bore 2AL?* Neg (a — By) 9" (a— By) (42) 
y 
a 
| / 
/ 
/ 
| | / 
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FIGuRgE 8.—Flow pest a curved wall. 


A STUDY OF SECOND-ORDER SUPERSONIC FLOW THEORY 


According to equations (38) and (36), the solution is 


JPN 


aa a olay eylg'(z—By)?+H(e—By) (48) 


‘Imposing the approximate second-order planar tangency 
condition (equation (22b)) gives 


(2—32EN 


Mace ae 


HOOO) (44) 


so that 


Hry=— }az)e"(z) +7 MXN —2) 


sar, eaae} aa) 





The complete second-order perturbation potential is therefore 


g= glzr—By)—e? \e(2— By)e'(x—By}+ 


Af aN {XN —2 
28 - yig (r — ~py)P+7 Gr? 3 


As 

iat (46) 
The same result can be found by solving equation (42) 
directly, using the impulse method (reference 15, p. 164). 

Flow quantities at the surface of the wall can be related 
to their values at the plane y=0 by expanding in Taylor 
series and discarding terms in ê. In this way the streamwise 
velocity perturbation at the wall is found to be given by 


H ; AEN 2 ar rena z 
pge e) g ee) (47) 
The pressure coefficient at the wall can now be calculated 
from equation (27) which, upon replacing N by its value from 


equation (11), gives 


)st* 


eg'(x Taoa E iA 


C= leet 48) 


ates 


This is the well-known result of Busemann (references 22 


and 23). To second order, the surface pressure coefficient 
depends only upon the local slope. 


ROLE OF CHARACTERISTICS 


It was pointed out previously that because of the under- 
lying significance of the characteristics for solutions of 
hyperbolic equations, it might be expected that the approxi- 
mate characteristics of the iteration equation would have to 
be revised successively at each stage. However, an iteration 
procedure was adopted which involves no such revision. 
It is therefore pertinent to mquire in this simple example 
What roles have been played by the actual and approximate 
characteristics. 

The flow past a single curved wall is given (until shock 
waves form) by a simple wave or distributed Prandtl-Mevyer 
expansion. Of the two families of characteristics, those of 
primary importance in a simple wave run downstream away 
from the wall. We therefore confine attention to that family. 

For the first-order equation, equation (32), the conven- 
tional theory (e. g., reference 15, ch. 5; reference 16, ch. 2) 
shows that the characteristics of the downstream family 
are the lines of slope 


(49) 
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This means that to first order the actual characteristics are 
approximated by the Mach lines of the undisturbed flow. 

For the second-order solution, a closer approximation to 
the characteristics could easily be found. It can be shown 
that if the first-order streamwise perturbation velocity at any 
point in the flow is u®, then the revised local values of 
Mach number and 6 are approximately ? 


Uos E 4-pN-1) oT | (50a) 
pea YM —1=6 14-64-15, | 6%) 


Combining this result with the first-order solution (equation 
(41)) shows that the revised Mach lines have the slope 


MEN ev) | 


it—3- é eg (x— 
However, because the iteration scheme adopted does not 
allow for such revision, these are not actually the character- 
istics of the second-order equation. Instead, the character- 
istics of equation (42) are still the original Mach lines of the 
undisturbed flow. 

Physically, the characteristics are lines slong which 
discontinuities in velocity derivatives are propagated, and 
this definition is completely equivalent to the mathematical 
one (reference 15, p. 297). Therefore, in the second-order 
solution given above, discontinuities in acceleration will 
occur along the original characteristics rather than, as they 
more properly should, along the revised characteristics. 

Suppose, however, that no such discontinuities occur. 
For flow past an isolated body the downstream character- 
istics are also lines along which the velocity is constant, 


(51) 


provided that shock waves do not appear. Setting 
d¢.= ourit + o,¢y=0 
(52) 
d¢,=oydr+o,dy=0 
it is seen that the velocity is constant if 
dy $u Ên (53) 


dr n hy 


For the second approximation (equation (46)) the velocity 
is constant along lines of slope 


dy _ ae N 
dx 2 





ce’ (26) | (54) 


which, according to equation (51), are the revised character- 
istics. Consequently, although the characteristics have not 
been revised in the mathematical sense, the solution behaves 
physically as if they had, so long as discontinuities do not 
occur. The question of discontinuities will be considered in 
the next section. 

The connection between the original and revised character- 
istics can be interpreted physically. The right-hand side of 
the Iteration equation may be regarded as representing the 

14s implied by the superscript, these are regarded as first-order values, because they are 


determined from the first-order solution. From this point of view, Af and § are zero-order 
quantities. 
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effect of a known distribution of supersonic sources through- 
out the flow field. The influence of this source distribution 
spreads downstream along both families of original character- 
istics. The resulting velocity changes are just such that to 
second order the velocitics become constant. along the 
revised rather than the original characteristics. 

Finally, it is interesting to note that the second-order 
potential is constant on lines which bisect the original and 
revised characteristics. For, setting 





do=¢.dz+¢,dy=0 (55) 
¢ is found to be constant along lines of slope 
ALN 
ou =z] 1 + ee eg e —sv)| (56) 


FLOW PAST A CORNER AND A PARABOLIC BEND 


A simple case in which discontinuities may occur is that of 
flow past a sharp corner. The exact solution is known to 
involye an oblique shock wave with attendant velocity dis- 
discontinuities for compression and a continuous Prandtl- 
Meyer fan for expansion. 

Denoting the tangent of the deflection angle by e, positive 
for compression (see fig. 4), the function g(x) appearing in 
equation (39) is 


0 «r<0 
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FIGURE 4.— Flow past a corner. 


From equation (46) the second-order perturbation potential 
is found to be. 


re 1 (58) 


o=—Z(e— By) + sa- By) — 
to the right of the line z=8y and zero to the left. Conse- 
quently, in either compression (e€>0) or expansion (e<0) 
the second-order potential suffers a discontinuous drop along 
the Mach line from the corner, of strength proportional to 
the distance from the corner. Such a discontinuity cannot 
be admitted, which indicates that the iteration process fails 
in this region. 

In the case of compression, the solution can be corrected 
by analytically continuing the perturbation potential up- 
stream until it can be joined continuously to the free-stream 
potential. (This is permissible since the line of discontinuity 
is not actually a characteristic.) From the result of equation 
(56) the juncture is seen to occur along the line from the 
corner which bisects the upstream and downstream Mach 
directions, as indicated in figure 5. The adjusted juncture 
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corresponds to a shock wave, for it is known that an oblique 
shock bisects the Mach lines to a first approximation (refer- 
ence 16, p. 354). In the case of expansion, this type of 
correction cannot be justified since it would involve con- 
tinuation of the free-stream potential across a true charac- 
teristic. Instead, a Prandtl-Meyer fan must be inserted. 


yh 


FIGURE 5.—Mach lines before and after adjustment of potential discontinuity. 


Evidently the iteration process is successful except within 
an angular region of order e lying near the Mach line from the 
corner. In particular, the pressure is given correctly every- 
where on the surface of the wall. 

It is enlightening to observe that the alternative method of 
iteration, in which the characteristics are successively re- 
vised, fails in the same region. The potential is double- 
valued over & fan-shaped region in the case of compression 
and is left undefined over a similar region in the case of ex- 
pansion (see fig. 6). The same artificial corrections are 
necessary to complete the solution. 


Fieure 6.—Second-order flow past a corner using revised characteristics. 


Consider next flow past a parabolic bend which is repre- 
sented by 


y=3 ex? = £20 _ (59) 


From equation (46) the second-order perturbation potential 
is found to be 





ALYN +1)—3 PN 
BO ae oyp MN aye ayy 


(60) 


The potential and also the velocities are continuous, so that 
the previous difficulties do not occur. The acceleration is 
discontinuous across the original characteristic z= By, which 
in this.case happens to be also a revised characteristic. 
However, a new complication arises. It is known that, in 
the exact solution for the compressive case, the character- 
istics form an envelope, as shown in figure 7. Inside the 
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cusp the potential is triple-valued (reference 16, p. 111), so 
that a shock wave must be inserted. This envelope must 
also arise in the second approximation since the character- 
istics are no longer parallel. However, the second-order 
potential given by equation (60) is single-valued, so that it 
cannot predict the formation of an envelope. Again the 
iteration process fails in a part of the flow field. 

It can be seen that the alternative iteration process, using 
revised characteristics, will produce an envelope. 


CONVERGENCE FOR PLANE FLOW 


The examples just considered indicate that the success of 
the iteration procedure should be carefully investigated. A 
step of an iteration process may be considered successful if, 
in some sense, it significantly improves the solution. In 


- ee ee 
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Fiter 7.—First- and second-order flow past a parabolic bend. 


particular, one Is interested in the success of the second- 
order solution. 

[t should be noted that a divergent process may be suc- 
cessful for many steps and that, on the other hand, con- 
vergence does not necessarily imply success. In practice, 
however, one would expect a convergent process to be suc- 
cessful. As used here, success is a subjective notion, not 
umenable to analysis. Consequently, only the convergence 
of the iteration procedure can be considered in any detai. 

Unfortunately, proofs of sufficient conditions for converg- 
ence have not been obtained, even in the case of plane flow. 
However, the preceding examples suggest certain conjectures 
regarding convergence. These will be stated and some argu- 
ments for their plausibility will be advanced. 

For flow past a slightly curved plane wall represented by 
y= (x), the solution obtained by iteration using the revised 
characteristics is conjectured to converge in any bounded 
region adjacent to the wall, provided that 

(a) e is sufficiently small, 

(b) gix) is continuously differentiable. 

If g(x) has only a piecewise continuous derivative, the 
convergence holds except possibly in fan-shaped regions 
springing from each corner, which lie near the original 
Mach line and subtend an angle of order e. 

For the iteration process actually adopted, in which the 
characteristics are not revised, the first n steps are conjec- 
tured to form part of a convergent process, provided that 

(a) eis sufficiently small, 

(b’) g(x) has continuous derivatives up to (n—1) st order 
if the potential is required, nth order if the velocities are 
required. 
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If condition (b’) is satisfied only piecewise, the result holds 
except possibly in fan-shaped regions springing from each 
Corner. 

In the first case, condition (a) is necessary in order to 
insure that the solution be unique, as is clear from the exam- 
ple of the parabolic wall. The preceding examples also 
show that condition (b} is necessary. i 

If the sufficiency of these two conditions is assumed, their 
connection with condition (b’) in the second case can be 
illustrated by analogy with a mathematical model which 
retains the essential difference between the two iteration 
processes; namely, that the correct characteristics are not 
used in the method actually adopted. Consider the first- 
order problem given by equation (40): 


Enr Pzr =À 
py (x, 0)}= eg (2) (61) 
P (0, y)=¢:(0, y)=0 


where we have taken 8=1 for convenience. The solution 


(equation (41)) is 

g=—egr—y) (62) 
Now we attempt to solve this problem using characteristics 
which differ from the true characteristics by Of{e). Thus 
we consider the equivalent problem 


Cy—(1 —e) Çr € rr 
Prlz, 0) = eg’ (x) (63) 
¢(0,y)=¢(0, y)=0 


and solve by iteration. In the first approximation the 
right-hand side is neglected, so that the differential equation 


becomes 
Pa? —(L—€) gx =0 (64) 
which has the solution, subject to the boundary conditions, 
p= —eg(r—vyl—ey) (65) 


Substituting this into the right-hand side of equation (63) 
gives the iteration equation for the second approximation: 


On —(1—€) ors =— êg" (z—y 1 —e€ 9) (66) 


Using the impulse method (reference 15, p. 164) gives the 
solution, subject to the boundary conditions, as 

ep = — eg(z—yl—ey) + 5 eyg (x—1—ey) (67) 
But this is just the Taylor series expansion, correct to 0(e), 
of the true solution (equation (62)). Subsequent iterations 
add additional terms to the Taylor series expansion. Hence, 
despite the use of slightly incorrect characteristics, the itera- 
tion process converges to the correct solution. The connec- 
tion between conditions (b) and (b’) is thus seen to be that 
the existence of sufficiently many continuous derivatives 
compensates for the fact that the wrong characteristics 
are used. 


i This model was suggested by Prof. C. R. DePrima of the California Institute of 
Technology. 
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AXIALLY SYMMETRIC FLOW 


Before discussing the general solution for bodies of revolu- 
tion, it is convenient to consider the simple problem of a 
cone. In this case the second-order solution can be found 
directly. The results will be useful in indicating which 
triple products should be retained in the general case. 


FLOW PAST A CONE 


Consider flow past a slender cone of semivertex angle 
tan! e, as shown in figure 8. The flow is conical and axially 
symmetric, so that the iteration equation is given by equation 





FIGURE &.—Flow past a cone. 


(17) with @ derivatives omitted. Including the boundary 
conditions from equations (19a) and (23a), the first-order 
problem. is 


a-a (0429-0 (68a) 
iade Abia (68b) 
` p=p =0 at i= œ (68c) 


The diferential equation can be solved by using the inte- 
grating factor t/y1— Ë. Imposing the upstream condition 
(equation (68c)) gives the well-known result 


g=—A (sech-! t— f1—&) . (69) 


which is understood to vanish except within the downstream 
Mach cone, whereiS1. The tangency condition (equation 
(68b)) is satisfied by putting 





2 


€ . 
AE e eh (Bo) (70) 
At the Mach cone (=1), all velocity perturbations vanish, 
so that the first-order solution predicts no deflection of the 
shock wave from the Mach cone (See reference 16, p. 403). 

Substituting the first-order solution into the iteration 
(equation (17)) gives 


sech! ¢ 1 
or F 





aip" 0+0- ae |2 (N—1) 


payiz] 


(71a) 
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and from equations (19b) and (23b) the corresponding 
boundary conditions are . 
at f= Be 


Bd’ = «(1 +¢—i¢’) (71b) 


(b=¢'=0 


Equation (71a) can again be solved using the integrating 
factor ¢t/¥1—#. The various integrals encountered can be 
treated by integrating by parts one or more times. Imposing 
the upstream condition (equation (71c)), the complete 
conical second-order perturbation potential is found to be 


at i= © (71c) 





$=—B (sech t— Y1—#)+A7M? | sech- a) — 


(NINTZE sech-! —4 l gra UZP tae tole (sech~ 5] 


(72) 


From equation (16), the streamwise and radial velocity 
perturbations are 


Y= —B sech-' t+ A*Af? | Goch f)?— 





U 
WOT 4 (N+)? ata VIZË z] (73a) 
pas gyt Pant] —22 a t 
OD HN 0 teaGtiyict BE =f) 


(73b) 


‘The constant B must be adjusted so as to satisfy the tangency 
condition (equation (71b)). In actual computation it is 
easier to. adjust B numerically in exactly this fashion rather 
than to calculate it from the cumbersome expression which 
could be written down. The pressure coefficient at any 
point can then be calculated from equation (26). 

The last term in the bracket in equation (71a) is the triple 
product Bpr pu which is retained in the iteration equation 
(equation (17)). Its retention is now justified by noting 
that its contribution—the last term in equation (72)—1s of 
the same order as other terms near the surface of the cone 
(¢=@e). Actually, it contributes a second term, which has 
been neglected since it is at most of order sech! fe. It 
can also be verified that the other triple products, whose 
form is indicated at the end of equation (17), are in tact 
negligible since they contribute at most terms of order 
é&(sech™! fe)’. Consideration of a further iteration indicates 
that a third approximation would add terms of order 
é(sech™! Be)’, which is greater than the terms just neglected. 

The second-order result for surface pressure coefficient is 
compared in figure 9 with the exact solution (reference 20) 
for cones of 5°, 10°, 15°, and 20° semivertex angles. Also 
shown for comparison are the first-order results based upon 
the exact isentropic expression (equation (26)) for the pres- 
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(b} Semi-vetex angle 10°. 
(d) Semi-vertex angle 20° 


Fiut BE §.—Comparlson of various Solutions for presure on a cone 


sure coefficient. The second-order solution is seen to provide 
a much better approximation throughout the range of 
Mach numbers up to the point at which the Mach angle 
equals the cone angle, beyond which the perturbation solu- 
tions have no physical meaning. 


SHOCK-WAVE POSITION FOR CONE 


The solution for plane flow fails near the Mach wave from 
a corner, Which suggests that the second-order solution for 
the cone may likewise fail near the Mach cone from the 
vertex. In the plane case, nevertheless, a first approxima- 


tion to the shock-wave position (and hence to the entropy 
change) can be calculated from the velocity perturbations 
near the Mach wave. We now consider whether this is 
true also for the cone. 

It was noted before that first-order theory predicts no 
disturbance at the Mach cone and hence no shock wave. 
According to second-order theory the velocity perturbations 
just behind the Mach cone are (equations (73)) 


~(0) a2 (0) 


—25f°Net (74) 
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so that the perturbation is normal to the Mach cone. Here 
“A (equation (70)) has been approximated by ¢&. From 
equation (50b) the cotangent of the revised Mach angle 
just behind the Mach cone is found to be 


BY = p[1—2AL‘N(N—-1) 4] (75) 


The upward stream inclination there is approximately 
(v/T7) 21, SO that the Mach lines have the slope 


dr_ 1 Í 
de (A+2MIN «cS . (76) 
Now if this can be taken to be the slope also of the revised 
characteristics just behind the shock wave, then the slope 
of the shock wave differs from that of the original Mach 
OEN ? M'N? 2 Me 
Í : 1 

tan oa ae ta = et 
This problem has been treated rigorously by Lighthill (refer- 
ences 24 and 25) and by Broderick (reference 26), who find 
that actually 


(77) 


_i_3 (y-+1)?MF t.. 
tan À a3 pe (78) 
which is one and one-half times the preceding result. The 


discrepancy means that the second-order solution fails near 
the Mach cone. The nature of this failure and the proper 
method of remedying it have recently been studied by Light- 
hill (reference 25). 

The entropy increase through a weak oblique shock wave 
is proportional to the cube of its inclination away from the 
Mach lines. Consequently, the entropy rise through the 
shock wave from a cone is O(c), as noted by Lighthill 
(reference 24). 


PARTICULAR INTEGRAL FOR AXIALLY SYMMETRIC FLOW 


Consider flow past a body of revolution which is either 
a slender pointed body with nose at the origin or one which 
extends indefinitely upstream with constant radius a for 
a0 (see fig. 10). The latter shape corresponds to the 
external flow past a sharp-edged, open-nosed body with 
supersonic flow at the lip. With slight modification the 
subsequent development can be applied to internal flow as 
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FIGURE 10.—Flow pest bodies of revolution. 


well. The meridian curve can be represented in the first 


case by 
r=R(r)=e9(z) " 220 (79a) 
and in the second by E 
a 230 
— = b 
PE T «20 Soa 


Here ¢ is again a parameter small compared with unity, and 
g(z) is a continuous function of order unity. 
The first-order problem is 


prt —Byr=0 (80a) 
p(z, R)}=R" (2) [1 + p:(2, W] (80b} 
(0, r)= p40, r)=0 (80c) 
The solution is known to be (reference 27) 
[F Aad 
an= eo 
cosh- Z>? = = aA l 

-f e F(z— gr cosh ujdu (81) 


The second form is useful for carrying out differentiation, 
after which the first form can be restored. The derivatives 
which will be required are 


cash-! ba i I—ßr 4 
a=- f fr "(2 —Br cosh ujdu= -f gA 
(822) 


I—b 
oah- -gr 
- e=8 | j f'(x—gpr cosh u) cosh u du= 
Q 


1 FA (a— FEE 
rje  (e—’— 67" 


In carrying out the differentiation the fact has been used 
that {(6)=0 for a body with finite slope. With coordinates 
as shown in figure 10, the lower limit of integration 6 is 0 
for the pointed body and — a for the semi-infinite body. 
The function f(z) may be interpreted physically (aside from 
a numerical factor) as the strength of a supersonic line source 
along the z axis. It is determined by the tangency condition 
(equation (80b)) which gives the following integral equation 
of the Volterra type for f": 


a (z-di on “teas 
0 (a — §)*— p72) 


R R' i £—~-f R(x) 
(z) ol -f JeF- BR) 


The second-order iteration equation is found from equa- 
tion (13) to be 


~ - (82b) 


| rr += — B’ pa = M?[2(N—1) Beever t+ 2010 tert 


Olor Prr, Pr Pr, zr Prr) (84) 
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The solution for the cone suggests that the terms indicated 
at the end of the equation are negligible. 

It will now be shown that a particular integral of this 
equation is given in terms of the first-order solution by 


| 1 
Krr)=M| ode+Nre)—zree | (85) 
The first group of terms contributes the first two terms on the 
right-hand side of equation (84), as can be verified by direct 


substitution. The last term in equation (85) accounts for 
the term gy, as follows: 


aetra ot) — gre?) = 
zet 12 Pr J5 


a 3 3 
Pr Err — Bor i rompra ren?) (86) 





-where repeated use has been made of the fact that ¢ satisfies 
equation (80a). The last group of terms consists of triple 
products involving z derivatives, which have already been 
neglected in equation (84), so that the result is proved. 

The complementary function x is a solution of equation 
(80) and can be written as 


oy (Pr Fede a E 
n=- | Que : F(x ie 


Using equations (82) the second-order tangency congitign 
{equation ee is found to be 


t-BR (x—BF"(OdéE 
ane coe 
f E—-BR FR 
Ry tin) | e (88) 


whieh is again & Volterra integral equation. 
METHODS OF SOLVING INTEGRAL EQUATION 


Discovery of a particular integral for bodies of revolution 
reduces the second-order problem to the same form as the 
first-order problem; namely, the solution of a Volterra in- 
tegral equation. Various methods of attacking this problem 
are listed by Hayes (reference 17, p. 140). 

An indirect method consists in assuming that the unknown 
source strengths in equations (81) and (87) can be represented 
by a few terms of a polynomial, for example, 


f2)=Cet+Oet+ ... +0," (89) 


The resulting solutions were introduced in a more formal 
manner by Haves (reference 17, p. 38), who has discussed 
their properties in detail. The first term alone gives the po- 
tential for the cone, (equation (69)). Additional terms give 
the solution for simple families of shapes. However, the 
method is not suitable for bodies having discontinuities in 
slope or curvature. Consequently, a more direct procedure 
is desirable. 

von Kármán first introduced an asymptotic solution of the 
integral equation (equation (83)) which has come to be 


known as the slender-body approximation (reference 27). 
For slender bodies, the source strength f(z) appearing in 
equation (81) is found to be approximately proportional to 
the rate of change of cross-sectional area. Thus 


fle) = 5=  [wRX2)] = R(2)B"(2) (90) 
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FIGURE 11.—Equivalence of polygonal source and sum of conical sources. 


Lighthill has shown (reference 19) that if R(x) and its first 
two derivatives are of order e and R’ is continuous, this de- 
termination of f(z) is correct to the order of terms retained 
in the first-order solution. For purposes of the second-order 
solution, it can be shown that f(x) may be determined in this 
way only if the first four derivatives of R are of order e and 
R” is continuous. This means that the body must have 
continuous curvature, which is a severe limitation. More- 
over, the slender-body approximation is found generally to 
cause unnecessary loss of accuracy even though the mathe- 
matical order estimate of the error is small. Consequently 
this approximation should be avoided if possible. 

The most satisfactory way of solving the integral equations 
is to use a step-by-step numerical procedure. In first-order 
theory the usual method, introduced by von Kármán and 
Moore (reference 28), is to assume that the unknown source 
distribution can be approximatel by a polygonal graph. 
This is equivalent to superimposing a number of conical 
source lines of different strengths, each shifted downstream 
with respect to its predecessor, as indicated in figure 11. 
The latter viewpoint is more convenient for computation. 
The strengths of the source lines are determined in succession 
by satisfying the tangency condition at a series of points on 
the surface of the body. The details of this prodedure are 
clearly explained in reference 3, page 77. 

For purposes of a second-order solution, this procedure 
must be modified in one respect. Unless the source distri- 
bution f(a) actually has corners, it must not be approximated 
by a polrgon. The reason is that a corner corresponds 
locally to adding a conical source line which would, according 
to the solution for the cone, produce false second-order dis- 
continuities in velocity and pressure across the Mach cone 
from the corner. Iristead, the procedure must be carried out 
using source lines of quadratic strength. The source strength 
f(x) can then be approximated smoothly so that false pressure 
jumps do not occur. <A single source line of this type repre- 
sents the flow past a slender pointed body with a cusped nose 
(see fig. 12), as is clear from the slender-body approxi- 
mation (equation (90)). s 
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METHOD OF SOLUTION FOR SMOOTH BODIES 


The second-order solution will be described first for bodies 
having continuous slope. Modifications for treating sharp 
corners will be discussed in the next section. 

The procedure is indicated in figure 13. The axis is 
divided into intervals by choosing points with abscissas 
Ex, at each of which a source line is to begin. Good accuracy 
is usually obtained if the interval length is not greater than 
8 times the local radius. The tangency condition will be 
imposed on the surface of the body at the points P,, which 
lie on the Mach lines from the points £y. 


—/ 


Ser 


FIGURE 12.— Body formed by source Hne of quadiatic strength. 
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FIGURE 18.—Method of solution for smooth bodies. 


For pointed bodies, the first-order solution is started with 
a conical source from the origin which gives the proper 
conical tip. This potential and the derivatives which are 
required are 


yo=— Cy (sech“t—y1—#) po => 
x 





=— Cy sech“! mE 
0, 0 Por T T= fi 
CEL (91) 
= BC, ~ , -o p= 


where 





Co Y1 — pett e? sech! Be 
Here Cy is the same as the A of equations (69) to (73), e 
being the tangent of the semivertex angle of the conical 
tip. No such term is required for the semi-infinite body. 
The subsequent procedure is the same for either body. 
Quadratic source lines are started from each of the points 
i, &, and so forth. For the pointed body & also lies at the 
origin, while for the semi-infinite body it is at —fa. For 


the nth such source line, the potential and its derivatives 
are given by 


Ya=—C,(z—§,F [G+ Ta : sech! ra —g v1 ar | 


2, =—2C, (x— Ey) —" — 1 — Fa) 


Pn, =60, (2—6) Vir a ae ci ra) 


Pa = — 2 Ca Sechi ry 





(92) 
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where 


The constants C, are determined successively by imposing 
the first-order tangency condition in turn at each of the 
points P,i;. From equation (80b), the condition is 

n n 

> Ox =R (i +>) Pn.) at Presi (93) 
where the summation begins with n=0 for the pointed body 
and n=1 for the semi-infinite body. In this way, values 
of the complete first-order potential » and its first and second 
derivatives are calculated at each of the points on the body. 

The velocities due to the particular second-order integral 

ý can then be calculated at the same points. Differentiating 
equation (85) gives 


VI | (e+Nro et edeetNron)—3 ro on | 


Yr =Af (p + N Ter) Per +¢,[(NV + lg, +N; Perr] ~ p lpr +3r¢e,;,) 
(94) 


Finally, the second-order complementary function x is 
determined by repeating the procedure used for ¢, finding 
new constants such that the second-order tangency condition 
is satisfied. According to equation (19b), the condition is 


VAS Sr iA) 


The second derivatives of x need not be calculated. 

The complete second-order perturbation velocities are 
found as the sums of the contributions from yẹ and x. Then 
the pressure coefficient can be calculated at each point P, 
from equation (26). 


at Pisi (95) 


TREATMENT OF BODIES WITH CORNERS 


Suppose that the meridian curve of the body has 2 sharp 
corner, which for convenience may be assumed to lic on the 
Mach cone from the origin, as indicated in figure 14. Then 
the method of solution must be modified for two reasons. 
In the first place, the intervals between source lines would 
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have to be chosen extremely small in order to obtain an 
accurate first-order solution bebind the corner. This diffi- 
eulty can be overcome by adding a new solution which causes 
a sharp deflection of the streamlines. In the second place. 
even if the first-order solution is determined exactly, the 
second-order solution does not yield the Busemann result 
just behind the corner, as it should since the flow is locally 


plane. This defect is remedied by properly canceling a 
r 
STET 





Figt'RE 14.— Body with corner. 


discontinuity which arises in the second-order solution at 
the corner. 

These two modifications require special solutions of the 
first-order equation which along the Mach cone from the 
origin have discontinuities in velocity in the first case, and a 
discontinuity in potential in the second case. Such solutions 
can be found by approximating to equation (83) in the vicin- 
ity of the Mach cone. Imposing the condition that there z 
is only slightly less than @# and keeping only leading terms 
in the difference r—8R leads to an Abel integral equation 
for the source strength. Inverting the integral equation 
shows that a potential having discontinuous nth derivatives 
results from a source distribution along the axis which is 


Setting f(z) =2"-7 in equation (81) 
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This mtegral represents the analytical continuation of the 
hypergeometric function, so that, except for a constant 


factor, 
r—fr 
“sir ) 


where a is the radius at the corner. The potential is under- 
stood to vanish except within the downstream Mach cone 
from the origin. The hypergeometric functions occurring 
here can all be expressed in terms of complete elliptic 
integrals with real moduli. 

In the first-order solution, e sharp deffection of the stream- 
lines at a corner is produced by adding a multiple of the 
potential which has discontinuous first derivatives. This is 
found by setting n=1, which gives 


g=(r— ary/* (p3 5 n+1;— (97) 
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Here t is the conical variable introduced in equation (14), 
and A and F are the complete elliptic integrals of the first 
and second kinds with modulus k=[(1—é)/(1+#)]”%. 
From the tangency condition (equation (80b)), it can be 
shown that in order to account for the corner the solution 
given by equation (98) should be multiplied by 


(By’— Ro) [1 Hiera] 
B+ ht,’ 


Here R,’ and Ry’ are the slopes of the meridian curve just 
ahead of and behind the corner and (¢.), is the value of or 
ahead of the corner. The first-order solution can thereafter 
be continued as described in the previous section. 

The second difiituity noted before was that the second- 
order solution is found to be incorrect just behind the corner. 
The proper method of treating this difficulty is to solve the 
case when the corner has been slightly rounded and then 
pass to the limit of a sharp corner. However, the following 
simpler procedure is found to give exactly the same result. 

The particular solution ý calculated from equation (85) 
is discontinuous along the Mach wave springing from the 
corner. If the discontinuity vanished at the corner, the 
solution could subsequently be revised as in the case of 
plane flow (see fig. 5). However, there is a finite jump in ¥ 
directly at the corner, which cannot be allowed. Conse- 
quently, the correction potential x must involve an equal 
and opposite jump. <A potential having such a discontinuity 
is obtained by setting n=O in equation (97). Then 


“ViVi | 
-N se E) 
-e EAG) 


Adding a suitable multiple of this potential cancels the dis- 
continuity in ¥. The second-order solution can then be 
continued as described in the preceding section. It can be 
verified that the pressure Jump at the corner has then the 
correct second-order value. 

It is instructive to analyze the behavior of a corner from 
another viewpoint. It was pointed out before that the right- 








(99) 


(100) 
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hand side of the iteration equation can be considered to 
represent the effects of a known distribution of sources 
throughout the flow field. In the case of a slightly rounded 
corner, this source distribution will be weak except between 
the Mach lines from the corner. As the corner shrinks to a 
point, the source intensity will increase in that region in 
such a way that the total strength remains constant. In 
the limit, the source distribution will behave like a Dirac 
delta function along the Mach line from the corner. The 
particular integral for plane flow (equation (35)) takes 
account of this impulsive function so that the correct solu- 
tion is automatically obtained. In the case of axially sym- 
metric flow, however, it is clear that the particular integral 
given. by equation (85) misses the contribution of the im- 
pulse. Itis therefore necessary to correct this shortcoming 
by adding the step-function potential given by equation 


(100). 
COMPARISON WITH NUMERICAL SOLUTIONS 


The accuracy of the second-order solution for bodies of 
revolution can be evaluated by comparison with examples 
calculated using the numerical method of characteristics. 

The first body to be considered is & circular-arc ogive of 
12%-caliber radius of curvature followed by a cylinder, 
which has a half angle of 16.26° at the tip. The second- 
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FIGURE 15.—Pressure distribution on smooth body of revolution at M=3.24. 


i ee ae ae we oe =" 


——— first order 
Second order E 
o Method of characterisiics 
(reference 3f)} 








Ficurs 16.—Pressure distribution on body of revolution with corner at Mm2.075. 


order solution has been calculated for this body at a Mach 
number of 3.24. This represents a severe test of the method 
because the Mach angle is then only 10 percent greater 
than the tip cone angle. Intervals were chosen such that 
the points P, lay at 0.1, 0.25, 0.5, 1, 2, and 3.5 calibers. 
The pressure distributions calculated by first- and seconid- 
order theory are compared in figure 15 with the results of 
various computations by the numerical method of character- 
istics. Of the latter, the result obtained from the interpola- 
tion chart given by Rossow in reference 29 is believed to be 
more accurate than the earlier German computations which 
were taken from the summary report of reference 30. Ex- 
cept near the tip, the second-order solution agrees very 
closely with the numerical results. 

The second body to be considered consists of a cone of 10° 
semi-vertex angle followed by a cylinder. The character- 
istics solution for this body at a Mach number of 2.075 hus 
been given by Liepmann and Lapin in reference 31. The 
first- and second-order solutions were calculated beyond the 
corner using the modifications discussed in the preceding 
section. Figure 16 shows the shape of the body, the loca- 
tion of source lines, and the pressure distributions caleulated 
by first-order theory, second-order theory, and the method 


-of characteristics. Again, the second-order results agree well 


with the characteristics solution. 


.. SERIES EXPANSION WITH RESPECT TO THICKNESS 


An alternative method of solving the exact perturbation 
equation (equation (4)) by successive approximations is to 
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ussume that the solution can be expanded in powers of the 
thickness parameter e Thus the exact perturbation poten- 
tial is written as 
(1) (2) (3} (4) 

$—Ge1621E6AiA+Gett _.. (101) 
Substituting into equation (4) and equating like powers of e 
yields a sequence of equations 
a) @) 
oe — 6 b2=0 
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which can be solved in succession. ‘The first is again the 
usual linearized equation. This method was applied to plane 
subsonic flow in references 6 and 12. 

Schmieden and Kawalki first pointed out (reference 14) 
that the power series assumed here does not always exist, 
even for plane flow. In general, terms of the form e*ln*e 
appear, beginning with elne in the third-order solution for 
plane flow and in the second-order solution for axially sym- 
metric flow. Furthermore, for a body of revolution only 
even powers of e arise. Consequently, it is necessary to 
assume a more general series of the form 


{1) {2) (3} (4) (5} 
p= pet pette in ettet... (103a) 


for plane flow and of the form 
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(103b) 
for axially symmetric flow. 

On the basis of this assumption, Broderick has developed 
a second-order solution for supersonic flow past slender 
pointed bodies of revolution (reference 32). The analysis 
is rather lengthy, since the simplification resultmmg from the 
discovery of a particular integral does not appear. The 
results are limited to shapes for which the cross-sectional 
area is given by an analytic function, or at any rate has its 
first four derivatives small of order «, and the first two con- 
tinuous. This is a severe limitation since, for example, 
the two bodies discussed in the previous section are not 
admissible. 

Broderick’s method vields the slender-body counterpart 
of the present second-order theory. Just as the usual first- 
order slender-body results can be derived as asymptotic 
forms of linearized solutions, so Broderick’s second-order 
slender-body results can be obtained by expanding the 
present second-order solution in powers of £ and łn t for 
small ¢ and retaining secondary as well as leading terms. 
The logarithmic térms arise from the series 

1 3 


sech i= Ina f— ., 
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(104) 
The expansion will now be carried out for the case of flow 
past a cone. 
It is clear from equation (70) that the constant A in the 
first-order solution (equation (69)) is given approximately by 


A=’... (105) 
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Substituting this value into equations (73a) and (73b), 


expanding in powers of ¢ and én ¢, and imposing the tangency 
condition (equation (71b)) shows that 


Bae+e| (2M—1)In MAN +1)—5 |+ -+ - (108) 


Then according to equation (73), the velocity perturbations . 
on the surface of the cone are 
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Replacing N by its value from equation (11), the approxi- 
mate pressure relation of equation (27) gives for the pressure 
coefficient on the surface of the cone 
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This is Broderick’s result (reference 32, equation (81)). 

This series is compared in figure 17 with the original form 
of the second-order solution which uses the exact pressure 
relation. For the most slender cone, the expansion in series 
causes only a moderate loss in accuracy. For more practical 
thicknesses, however, the expansion reduces the accuracy to 
such an extent that for the cone of 20° semivertex angle 
Broderick’s solution is inferior to the first-order result (with 
the exact pressure relation). The reason must be that the 
iteration process itself converges more rapidly then do the 
subsequent expansions which are required to reduce it to the 
slender-body series form. Hence terminating all expansions 
at terms of the order of those retained im the iteration process 
results in an unnecessary loss of accuracy. 


(108) 


THREE-DIMENSIONAL FLOW 
PARTIAL PARTICULAR INTEGRAL 


It might be hoped that a particular integral, which so 
greatly simplifies the iteration for plane and axially sym- 
metric flows, could be found for the general three-dimensional - 
case. The various methods of existing first-order theory 
could then be applied immediately to the problems of second- 
order flow past such shapes as inclined bodies of revolution 
and three-dimensional wings. 

A part of such a particular integral is found at once, being 
common to the two special cases. Consider the three- 
dimensional iteration equation (equation (12)), which may be 
written 
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Exact (reference 20) 
Second order l i 
Second order exponded in series. 
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FIGURE 17.—Eflect of expanding in series upon second-order pressure on a cone, 


It can be readily verified that except for the term in N and 
the triple products (the last three terms) a particular integral 


Is given by 
Y =M topz 
which appears in both equations (35) and (85). 
The iteration equation is thereby reduced to 


Pr -+ Pr:— P brs =M? (2 PN PrPrz + Pr Pvy +2 PyLl:P re + Pr Psr) 
(111) 


(110) 


It hus not been possible to find a particular integral of this 
equation in terms of the first-order potential. The solutions 
for plane and axially symmetric flow do not appear to suggest 
a generalization. On the other hand, there is no assurance 
that such an integral cannot be found. When the triple 
products are negligible, the right-hand side of equation (111) 
vanishes for y=—1 (N=0). However, investigation of the 
previous solutions indicates that the idea of here taking 
y=—! is not legitimate. 
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In the absence of a complete particular integral, the re- 
duced iteration equation (equation (111)} must be attacked 
by more conventional methods. In principle, it is always 
possible to find a particular integral of a linear nonhomoge- 
neous equation with the aid of the fundamental solution 
associated with the differential operator. For the tbree- 
dimensional wave operator which occurs here, the funda- 
mental solution is 


—Biy—aPt+le— t] 
which can be interpreted as the potential at any point 
(z, y, z) lying inside the downstream Mach cone from & unit 
supersonic source at ($, n, ©). With the aid of Green’s 
formula, it can be shown that a particular integral of 


E 1 
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Prr FOr B brs = F(z, Y, 2) (113) 
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where the integration extends throughout that portion of 
the forward Mach cone from the point (x, y, z) within which 
F is defined. 

In practice, the integration indicated in equation (114) is 
generally not feasible. For example, even the simplification 
of axial symmetry reduces equation (114) only to a double 
integral of F(z,r) multiplied by a complete elliptic integral 
of complicated argument. Avoiding such integrals by dis- 
covery of the particular solution clearly represents a great 
simplification in this case.$ 

In the following sections, one example of a three-dimen- 
sional solution will be given, and the possibility of treating 
other shapes will be discussed thereafter. 


FLOW PAST AN INCLINED CONE 


The problem of a cone at an angle of attack illustrates the 
use of separation of variables to reduce the three-dimensional 
iteration equation to tractable form. 

Two alternative coordinate systems are suitable for bodies 
of revolution at an angle of attack. In wind axes the body 
is inclined, while in body axes the stream impinges en the 
body obliquely. The latter system is simpler for first-order 
problems and is probably better for the second-order solution 
also. However, wind axes will be used here, since otherwise 
the iteration equations must be rederived. 

To simplify the solution, it will be assumed that the angle 
of attack æ is so small that its square can be neglected. This 
will give a solution nonlinear in the body thickness but linear 
in æ, and will therefore yield the initial slope of the lift curve 
correct to second order. The coordinate system is indicated 
in figure 18. 


i Comparing the two methods would lead to the evaluation of definite Integrals Involving 
complete elliptic integrals, which might be of some interest. 

t In previous versions of this work (references I and 2) an oblique transformation was ap- 
plied which effectively unyaws the axis of the body. However, this additional transforma- 
tion has been found to complicate rather than simplify the analysis, and has been omitted 
here. 
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To this approximation the surface of the cone is given by 
r=R(r, 6)=[e—a(1+e”) cos Az 
i= T(6) =B [e—a(1 +e) cos 8] 


(115a) 
(115b) 


The first-order equation for the conical perturbation poten- 
tial is given by the left-hand side of equation (17): 


(11. +s =0 (116) 


The solution required here is the sum of potentials for a coni- 
cal line source and dipole (reference 3, p. 74) and has the 


form 
o=—A( sech-¥—+f/1—-F a cos @ 
(117) 
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From equations (16) and (19a) the tangency condition is 
found to be 


Bo:=[e—a(1+e*} cos 6](1+¢—te,) at t=T (118) 


Substituting from equation (117) and expressing values on 
the cone in terms of their values at i= fe by means of Taylor 
expansions, it is found that 


FIGURE 18.—Inclined cone In wind axes. 
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According to equation (110), a partial particular integral 
of the second-order iteration equation is, in conical form, 


pa = MGi d= AMH A sech t(sech— (—./1—#) — 
—— f= -1 
o| 311P sech 1—# sech t_,i—t a E E 0? |e ae = 


Pot tia cos 8 (120) 


There remains to solve the reduced iteration equation given 
by equation (111), which becomes 
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Exact (reference 35) 
——-—— Second order expanded in series 
-~-----~ First order {Tsien} 
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FIGURE 19.—Comparison of various solutions for normal force on a cone. 


This is reduced to two total differential equations by setting 


H(t, == do(t)+¢,(t)a cos 8 (122) 


Here ¢ is associated with the axial component of the free 
stream, and $, with the cross flow. The first of these is 
known from the previous solution for the symmetrical cone 
(equation (72)). The equation for the cross-flow term ¢, is 


n$ 2 sech! 
Op mA OM (w at 
NFPA vi) (123) 
Setting 
$ (t)=te (t) (124) 


reduces this to a linear first-order equation in w which can be 
integrated to find that 
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The tangency condition separates into the two conditions 


Bo! + $0) ef + otp O] ati=8e (126a) 
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The first of these is the previous relation (equation (71b)) 
which determined the constant B in equation (72). Simi- 
larly, the second of these determines the constant D in 
equation (125). 


4 
SERIES EXPANSIONS FOR PRESSURE AND NORMAL FORCE ON INCLINED 
Numerical results have been calculated only for the case 
in which the solution is expanded in powers of ¢ and in ¢ (the 


slender-body approximation). Carrying out the expansion, 
the constant D is found to be 


D=26¢ E +2Mtein j (3MN+5) ae... ‘| (127) 


Then caleulating the velocity components from equation 
(117) and the pressure from equation (27) gives, on the 
surface of the cone, 


CO ke E er (Min a3 M+1) | oh ae 


(128) 
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Here (Cy)o is the value for zero angle of attack, given by 
equation (108). Integrating gives the normal-foree coeffi- 
cient, based on cross-sectional area: 


Cyne e tafe (Atn 5-5 Me+1)+ hs a 
5 poU*(area) j 
(129) 


This result has been obtained also by Lighthill (reference 33), 
who has calculated the lift on bodies of revolution having 
analytic meridian curves by assuming a series expansion for 
the velocity potential. 

Stone (reference 34) has developed a solution for inclined 
cones Which is linearized with respect to ea, but otherwise 
exact. Kopal (reference 35) has published tables of the 
numerical results of Stone’s theory. A comparison of eque- 
tion (129) with this exact theory and with Tsien’s first-order 
solution (reference 36) is shown in figure 19 for 5° and 10° 
cones. The earlier discussion of series expansions suggests 
that the agreement might improve if the solution were not 
expanded in series. 


SHOCE-WAVE POSITION FOR INCLINED CONE? 


Just behind the Mach cone the velocity components are 


uy _ife —9 IN 4 | = 
-( p) a(t), 222. et(L—88æ cos 6) (130) 
For simplicity, using equations (119), A and C have here 
been approximated by eand 28e?. Comparing equations (74) 
und (77), it is seen that if these were the velocities just 
behind the shock wave, then the difference between the 
shock-wave angle and the [ach angle would be 
h—sin! Ee MN (18a cos 8) (131) 
Hence the ratio of the angular rotation of the shock wave 
to that of the cone would be 
- = §B2VL2N%et (132a) 
Recently Lighthill has derived a simple expression for the 
shock-wave position for any conical body lying inside 
the Mach cone (reference 25). From his results it is found 
that the preceding relation is incorrect, the correct expression 
being again one and one-half times as large, so that 
§ AL , rs 
AT 1)? ar e*== 12 M BIN? et (132b) 
For moderate cone angles this latter result agrees well with 
the exact values calculated by Kopal (reference 35) from 


Stone’s theory. Again the discrepancy indicates failure of the 
second-order solution near the Mach cone. 


EXTENSIONS OF THE THEORY 


Two important classes of problems have only been touched 
upon here. One is wings; the other, inclined bodies of 
revolution. The problem of inclined bodies has recently 
These relations were given incorrectly in references 1 and 2. 
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been studied further in reference 37. The iteration equations 
were there rederived in body rather than wind axes, which 
simplifies the tangency condition. Again only a partial 
particular integral could be found. 

The possibility of discovering particular integrals of the 
iteration equation might be investigated more systematically. 
If none can be found for general three-dimensional flow, 
special cases such as conical flow should be studied. 


TREATMENT OF WINGS 


One of the most useful applications of first-order theory is 
to thin flat wings. For conical wings, the reduced iteration 
problem can be transformed, by the standard conical theory 
(reference 18}, into the problem of solving Poisson’s equation 
inside a circle. This case has recently been studied by 
Moore (reference 38), who has calculated results for the 
nonlifting wing of triangular plan form lying inside the Mach 
cone and of symmetrical wedge section. Although the two 
extremes of the plane airfoil and the slender cone show that 
compressive pressures are reinforced in the second approxi- 
mation, Moore calculates a reduction for an intermediate 
sweptback case. 

Two difficulties can be anticipated. First, if the wing has 
subsonic edges, infinite velocities arise there, so that the 
assumption of small perturbations is violated. It is known 
that in first-order theory this is no essential objection, 
since the pressure is found correctly except in the immediate 
neighborhood of the singularity, and the integrated values of 
lift and moment are correct to first order. Schmieden and 
Kawalki (reference 14) and Kaplan (reference 13) have 
indicated that this result extends to the second approximation 
for subsonic flow, so that probably no real difficulty exists. 

Secondly, if the wing has supersonic edges, the failure of 
the iteration process along Mach lines from the apex can be 
expected to affect the surface pressures. Again it is possible 
that integrated values will be correct to second order. 
Otherwise, it may be possible to adjust the solution in those 
regions, in a Manner similar to that shown in figure 5. 


HIGHER APPROXIMATIONS 


It seems unlikely that third or higher approximation 
would ever be justified. Other neglected factors, chiefly 
viscosity and heat conduction, should ordinarily be considered 
first. However, the Busemann second-order result has been 
extended to third and even fourth order (reference 39), 
and various writers have considered the third approximation 
for plane subsonic flow (references 7,9, and 11). Ifa third 
approximation should be considered worthwhile, the iteration 
could be repeated. Again the cases of flow past a curved 
wall and a cone would serve as helpful examples. 


APPLICATION TO SUBSONIC FLOW 


The iteration equation and the particular integrals are in 
no way restricted to supersonic flow. The particular solution 
for plane flow might profitably be compared with the sub- 
sonic solutions of references 5 to 4. 

Recently Kaplan has shown (reference 40) that the par- 
ticular integral for plane flow (equation (35)) can be derived 
formally, using complex variable theory, and has similarly 
obtained the third-order particular integral. The latter is 
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not applicable to supersonic flow except in special cases which 
are free of shock waves. 

The particular integral for plane flow has recently been 
applied to several subsonic problems by Harder and Klunker 
(references 41 and 42). 

The particular solution for axially symmetric flow makes 
possible a second-order solution for bodies of revolution at 


subsonic speed. In this case, the integral equation can be 
treated by the methods used for the airship problem. 


AMES AERONAUTICAL LABORATORY 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Morrett Fisitp, CALIFORNIA, April 2, 1952 


APPENDIX 
SYMBOLS Ex abscissa of origm of nth source line 
a constant reference radius for body of revolu- į P local density E 
tion Ta conical variable referred to 2=&, as origin 
b abscissa at which source distribution for body rather than z=0 . l 
of revolution begins 9 first-order (linearized) perturbation potential, 
A, B,C,D constants determined by boundary conditions same as Ẹ™ 
c local speed of sound d second-order perturbation potential, same 
Ox constant coefficients of series as p” 
Cy normal-force coefficient p exact pertubation potential 
Cs pressure coefficient a nth-order perturbation potential 
E complete elliptic integral of the second kind | $ nth term in series expansion of perturbation 
with modulus k=[(1—2)/(1+0)]? ` potential | 
F(z), F(a) source-distribution functions for body of | x complementary function for second-order iter- 
revolution ation equation 
Falz, Y, 2) known right-hand side of (n+1) st-order iter- | ¥ particular integral of second-order iteration 
ation equation equation l 
gl) continuous function of order unity which | Ys partial particular solution for three-dimen- 
vanishes for z=0 sional flow 
hj, H, J arbitrary functions of one variable (See equation (110).) . 
complete elliptic integral of the first kind with | 28 complete velocity potential 
modulus s=[(1—#)/(1+4)]'? a auxiliary variable 
(See equation (124).) 
AM free-stream. Mach number (=) 
SUPERSCRIPTS 
2 
N ie oe = = (n) result of nth-order solution 
2 BF B ; f 
, — value in conical form 
p local statie pressure : (See equation (15).) 
Pi points on body of revolution at which tau- 
gency condition is imposed Sere 
g local speed of flow 0 associated with axial flow 
r radius in cylindrical coordinates 1 associated with cross flow 
R radius of meridian curve of body of revolution | ® ‘ free-stream conditions 
S function defining surface of body a value just ahead of corner 
bl Br b value just behind corner 
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